THE CONVEX HULL OF A VARIETY 



KRISTIAN RANESTAD AND BERND STURMFELS 

Abstract. We present a characterization, in terms of projective biduality, for the hyper- 
surfaces appearing in the boundary of the convex huh of a compact real algebraic variety. 



1. Formula for the Algebraic Boundary 

Convex algebraic geometry is concerned with the algebraic study of convex sets that arise 
in polynomial optimization. One topic of recent interest is the convex hull conv(C) of a 
compact algebraic curve C in R". Various authors have studied semi definite representations 
[HIIIH], facial structure [171 [23], and volume estimates [3l[T9] for such convex bodies. In [Tl] 
we characterized the boundary geometry of conv(C) when n = 3. The boundary is formed 
by the edge surface and the tritangent planes, the degrees of which we computed in [HI 
Theorem 2.1]. Here, we extend our approach to varieties of any dimension in any M". 

Throughout this paper, we let X denote a compact algebraic variety in which affinely 
spans M". We write X for the Zariski closure of X in complex projective space CP". Later we 
may add further hypotheses on X, e.g., that the complex variety X be smooth or irreducible. 

The convex hull P = conv(X) of X is an n-dimensional compact convex semialgebraic 
subset of M". We are interested in the boundary dP of P. Basic results in convexity [T0| 
Chapter 5] and real algebraic geometry P'j Section 2.8] ensure that dP is a semialgebraic 
set of pure dimension n — 1. The singularity structure of this boundary has been studied 
by S.D. Sedykh [201 Our object of interest is the algebraic boundary daP, by which we 
mean the Zariski closure of dP in CP". Thus daP is a closed subvariety in CP" of pure 
dimension n — 1. We represent daP by the polynomial in M.[xi, . . . , Xn] that vanishes on dP. 
This polynomial is unique up to a multiplicative constant as we require it to be squarefree. 
Our ultimate goal is to compute the polynomial representing the algebraic boundary daP- 

We write X* for the projectively dual variety to X. The dual variety X* lives in the 
dual projective space (CP")^. It is the Zariski closure of the set of all hyperplanes that are 
tangent to X at a regular point. According to the Biduality Theorem of projective geometry, 
we have (X*)* = X. We refer to [7i §1.1.3] for a proof of this important result. 

For any positive integer k we let X'*^! denote the Zariski closure in (CP")"^ of the set of 
all hyperplanes that are tangent to X at regular points that span a (fc— l)-plane. Thus 
X'^l = X* is the dual variety. We consider the following nested chain of algebraic varieties: 

xw c ■ ■ ■ c x[2i c xw C (CP")^. 
Our objects of interest is the dual variety, back in CP", to any X''^! appearing in this chain. 
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To avoid anomalies, we make the assumption that only finitely many hyperplanes are 
tangent to X at infinitely many points. Without this assumption, the relevant projective 
geometry is much more subtle, as seen in the recent work of Abuaf [Ij. With this assumption, 
for small values of k, the dual variety (Xt'^')* equals the k-th secant variety of X, the closure 
of the union of all {k — l)-dimensional linear spaces that intersect X in at least k points. 

The codimension of this secant variety is at least 2 ii k < L dim(x)+i J • ''^i-^) the 
minimal integer k such that the k-th secant variety of X has dimension at least n — 1. 
Thus we have r{X) > [ ^^^(xj+i l • "^^^ inequality k > r(X) is necessary for (Xt^^)* to be a 
hypersurface. The main result in this article is the following formula for the convex hull. 

Theorem 1.1. Let X be a smooth and compact real algebraic variety that affinely spans 
M", and assume that only finitely many hyperplanes in CP" are tangent to the corresponding 
projective variety X at infinitely many points. The algebraic boundary of its convex hull, 
P = conv(X), is computed by biduality as follows: 

n 

(1.1) daP C U (XW)*. 

k=r{X) 

In particular, every irreducible component of daP is a component o/(Xt^l)* for some k. 

Since daP is a hypersurface, at least one of the (X'^^)* must be a hypersurface. However, 
others may have higher codimension and these can be removed from the union. The reason 
for possibly not having equality in f ll.ip is that some hypersurface component in (X'^^l)* may 
have no real points, or its real points may be disjoint from the boundary of P = conv(X). 
Such components must also be removed when we compute the algebraic boundary daP- 

When the inclusion X''^! C X^^~^'^ is proper, the former is part of the singular locus 
of the latter. In particular X^^l is in general part of the /c-tuple locus of the dual variety 
Xm = X*. However, the singular locus of X* will have further components. For example, 
the dual variety of a curve or surface in CP^ has a cuspidal edge defined, respectively, by the 
osculating planes to the curve, and by planes that intersect the surface in a cuspidal curve. 

Our presentation is organized as follows. In Section 2 we discuss a range of examples 
which illustrate the formula fll.ip . The proof of Theorem 11.11 is given in Section 3. We also 
examine the case when X is not smooth, and we extend Theorem 1 1.1 1 to varieties whose real 
singularities are isolated. Section 4 features additional examples. These highlight the need 
to develop better symbolic and numerical tools for evaluating the right hand side of fll.ip . 

2. First Examples 

2.1. Polytopes. Our first example is the case of finite varieties, when dim(X) = 0. Here 
P = conv(X) is a full- dimensional convex polytope in M", and its algebraic boundary daP 
is the Zariski closure of the union of all facets of P. The formula (11.10 specializes to 

daP C (X["l)*. 



THE CONVEX HULL OF A VARIETY 



3 



Indeed, X'"! C (CP")^ is the finite set of liyperplanes tfiat are spanned by n affinely in- 
dependent points in X. Typically, this includes hyperplanes that do not support dP, and 
these should be removed when passing from (X^"!)* to daP- It is important to note that the 
Zariski closure, used in our definition of the algebraic boundary daP, depends on the field 

C M we are working over. If we take ii" = M then daP is precisely the union of the facet 
hyperplanes of P. However, if K is the field of definition of X, say K = Q, then daP usually 
contains additional hyperplanes that are Galois conjugate to the facet hyperplanes. 

Here is a tiny example that illustrates this arithmetic subtlety. Let n = 1 and take X 
to be the variety of the univariate polynomial — 3x + 1. This polynomial is irreducible 
over Q and has three real roots. The smallest root is a = —1.3888... and the largest root is 
P = 1.2146.... Clearly, P = conv(X) is the line segment [a, /3] in R^. If we take K = M. then 
daP = {a, /3}, but if we take K = Q then daP consists of all five complex roots of f{x). 

2.2. Irreducible Curves. Let n = 2 and X an irreducible compact curve in of degree 
d > 2. Since X is a hypersurface, we have r{X) = 1. Suppose that the curve X does 
not bound a convex region in R^. The algebraic boundary of the convex set P = conv(X) 
consists of X and the union of all bitangent lines of X. In symbols, 

daP C (XW)* U (xt^'i)* = X u (x['i)*. 

For a smooth curve X of degree d, the classical Pliicker formulas imply that the number of 
(complex) bitangent lines equals {d — 3){d — 2)d{d + 3)/2. Hence, daP is a curve of degree 

{d-3){d-2)d{d + 3) 
deg{daP) < d + . 

Next consider the case where n = 3, dim(X) = 1, and r(X) =2. If X is irreducible then 

daP C (X^l)* u (X[3])*. 

The first piece (X'^^)* is the edge surface of X, and the second piece (X^'^l)* is the union of 
all tritangent planes. For a detailed study of this situation, including pretty pictures of P, 
and a derivation of degree formulas for (X^^l)* and (Xl^l)*, we refer to our earlier paper [14]. 
Further examples of space curves are found in Subsection 4.1 below and in [TU Section 4]. 

Sedykh and Shapiro [19] studied convex curves X C R" where n = 2r is even. Such a curve 
has the property that |X fl if| < n for every real hyperplane H. The algebraic boundary of 
a convex curve is the hypersurface of all secant (r — l)-planes. In symbols, daP = (X'*"!)*. 

2.3. Surfaces in 3-Space. Let X be a general smooth compact surface of degree d in 
R'^. Confirming classical derivations by Cayley, Salmon and Zeuthen [16j p. 313-320], work 
on enumerative geometry in the 1970s by Piene fl3{ p. 231] and Vainsencher [221 P-414] 
establishes the following formulas for the degree of the curve X'^^, its dual surface (X^^l)*, 
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and the finite set Xl^l in (CP^)^: 



deg(X 



d{d -l){d- 2)(c/3 -d'^ + d- 12) 



2 

deg((X[2])*) = did-2){d-3){d^ + 2d-A), 

deg(X[=^]) = deg((X[3l)*) 

_ d^ - 6d^ + 15rf^ - 59d^ + 204:d^ - 339rf^ + 770^^ - 2056^^ + 1920c/ 

6 ■ 
We can expect the degree of daP to be bounded above by d plus the sum of the last two 
expressions, since 

daP c (xW)* u (x[2])* u {x^^^y = X U (X[21)* u (x[=^])*, 

unless X is convex or otherwise special. For a numerical example consider the case d = 4, 
where we take X to be a compact but non-convex smooth quartic surface in M^. The above 
formulas reveal that the degree of the algebraic boundary daP can be as large as 

deg(X) + deg((X[2])*) + deg((X[3l)*) = 4 + 160 + 3200 = 3364. 

2.4. Barvinok-Novik curve. We examine the first non-trivial instance of the family of 
Barvinok-Novik curves studied in [01 123] ■ This is the curve X C parametrically given by 

(ci, C3, si, S3) = (cos(6'), cos(3^), sin(^), sin(3^)) . 

We change to complex coordinates by setting xj = cj + ■ Sj and Xj = Cj — 1 ■ Sj. The 
convex body P = conv(X) is the projection of the 6-dimensional Hermitian spectrahedron 

/ 1 Xi X2 Xs\ 

Xi 1 Xi X2 
X2 Xi 1 Xi 
\X3 X2 Xi 1 J 

under the linear map MP — )■ R'*, (ci, C2, C3, si, S2, S3) (ci, C3, si, S3). The curve X is the 
projection of the curve in that consists of the above Toeplitz matrices that have rank 1. 

The convex body P = conv(X) was studied in pTf Example 5.5]. It is the 4-dimensional 
representative of the Barvinok-Novik orbitopes (cf. [6l[23]). Its algebraic boundary equals 

daP = (Xt'l)* U (X[31)*. 

The threefold (X'^^)* represents the 2-dimensional family of edges of P, while the threefold 
(X'^1)* represents the 1-dimensional family of triangles in dP, both of which are described 
in [HI Thm. 4.1]; see also |23]. The defining polynomials of these two hypersurfaces in Mf^ are 

( Xgxf — 2xfx3xfx3 + x^xl + Axfxl — 6xiX3xf — 6xfxiX3 + 12x^x3x^x3 from fX^^' V 

— 2x|xfx3 — 2x^X3X3 — 3x1x1 + 4x3xf + 4xfx3 — 6X1X3X1X3 + X3XI ) 



,Ci,C2,C3,Si,S2,S3) G 



is positive semidefinite > 
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and {X3X3-I) = (cl + sl-l) from {X^^^)*. 

Both of these threefolds are irreducible components of the ramification locus that arises when 
we project the hypersurface of singular Toeplitz matrices from into as above. 

3. Proof of the Formula 
We turn to the proof of our biduality formula for the algebraic boundary of P = conv(X). 

Proof of Theorem We first prove that the supporting hyperplane of any exposed face F 
of P lies in X^^^ for some k. Suppose that dim(F) = k — 1 and let Lp he the projective span 
of F. By Caratheodory's Theorem, every point of F lies in the convex hull of k distinct 
points on X. In particular, the {k — l)-plane Lp intersects X in at least k points that span 
a (fc — l)-simplex in F. If if is a supporting hyperplane for F, then H contains F and is the 
boundary of a halfspace that contains X. Since X is smooth, the tangent plane to X at each 
point qEXnFCXnH must therefore be contained in H. We conclude that [H] G X^''\ 

Now, consider any irreducible hypersurface Y C CP" whose real locus has full-dimensional 
intersection with the boundary dP C M". We need to show that y is a component of (X'*^])* 
for some k. In the next paragraph we give an overview of the proof that follows thereafter. 

First, we shall identify the relevant number A; = fcy + 1, by the property that Y has a 
linear space of dimension ky through every point. In fact, we shall prove that Y contains a 
unique fcy-plane through a general point of Y. Thus, at a general point, the hypersurface 

Y is locally a fibration. In particular, the general point in dP fl Y lies in a /cy-plane that 
intersects P along a fcy- dimensional face. Subsequently, we will show that the supporting 
hyperplanes of these faces are tangent to Y along these /cy-planes, before we prove that 
Y* C X'^l. From this, we shall finally conclude that F is a component of (X^'^^)*. 

Let g be a general smooth point in the (n — l)-dimensional semialgebraic set dPCiY. Since 
the union of the exposed faces of P is dense in dP, there exists an exposed face Fg that has 
q in its relative interior. The hypersurface Y contains the boundary of P locally at q, and 
hence it contains the face Fg. Since F is a variety, it contains the projective span Lp^ of the 
face Fq. Let fcy = dim{Lp^). Since g is a general smooth point in dP fl Y, the hypersurface 

Y contains a /cy-plane through every point of Y. In fact, since Fq is an exposed face, it is 
the unique face through q, so Y contains a unique fcy-plane through every general point. 

Next, let if be a hyperplane that exposes the fcy-dimensional face Fg of P. We will show 
that H coincides with the tangent hyperplane Hg to Y at q. As g is a general interior point 
in Fg, we then conclude that H is tangent to Y along the entire fcy-plane Lp^. The key to 
our argument is that H is assumed to be tangent to X at the points X Ci Fg that span Lp . 

If Y = Lp^ is itself a hyperplane, there is nothing to prove, except to note that fcy = n — 1, 
that H = Hg, and that Y* is an isolated point in Xt^^l Otherwise, we compare H and the 
tangent plane Hg via a local parameterization of Y at q. Let k = ky + ^ and m = dim(X), 
let pi, ..■,Pk be points in X fl that affinely span Fg, and let 



(for i = 1, ...,k) 
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be local parameterizations of X at the points pi. The point q lies in the affine-linear span 
of the points Pi, so q = Yl^ ^^iPi for some real coefficients Oj with = 1. There may 

be polynomial relations in the local parameters ti defining /c-tuples of points in X whose 
affine-linear span lies in Y. These relations define a subvariety Z in the Cartesian product 
X'' that contains the /c-tuple (pi, . . . ,Pk)- A local parametrization of Lp^ at g has the form 

with affine-linear functions ctj in the Ui. Since Y is locally a fibration, the algebraic functions 
7i and a provide a local parameterization of the complex variety Y at the point q: 

{u,ti,...,tk) ^ a{u) + Y^^^^ai{'yi{ti)) + e{u,ti,...,tk) 

Here, the function e only contains terms of order at least two in the parameters. The tangent 
space Hq at q is spanned by the linear terms in the above parameterization. But these linear 
terms lie in the span of {ai{u), q;„(m)) and the linear terms in (71, ...,jk)- The former span 
Lp^ , while the latter span the tangent spaces to X at each of the points p,. So, by assumption 
they all lie in the hyperplane H that supports dP at Fq. Therefore, the hyperplanc H must 
coincide with the tangent plane Hq to Y at q. Since g is a general point not just in Y but 
also in Lp^, we conclude that H is tangent to Y along the entire plane Lp^. 

We have shown that the tangent hyperplanes to Y are constant along the /cy-planes con- 
tained in Y. This imphes that the dimension of the dual variety Y* is equal to n — k where 
k — ky + l- Locally around the point q, these tangent hyperplanes support faces of dimension 
ky = k — 1 the convex body P. This ensures that the inclusion Y* C X^'^'^ holds. 

We next claim that Y* is in fact an irreducible component of the variety X'^l This will 
be a consequence of the following general fact which we record as a lemma. 

Lemma 3.1. Every irreducible component o/X^^l has dimension at most n — k. 

Proof. Let W C X^^l be a component, and let kw be the minimal / such that W is not 
contained in X''"*"^'. Then kw > k and is a component of X^''^^. 

Let CX C CP" X (CP")^ be the conormal variety of X, the closure of the set of pairs 
{p, [H]) e CP" X (CP")^ such that the hyperplane H is tangent at the smooth point p e X. 
It has dimension n — 1. By assumption, the projection p : CX (CP")^ into the dual space 
has only finitely many infinite fibers. Therefore X* = p{CX) is a hypersurface and W is 
part of its fc^y-tuple locus. If [H] is a general point in VF, then X* has at least kw branches 
at [H]. Let (pi, [H]), {pk^, [H]) be smooth points in CX in the fiber over [H], such that 
Pi,---,Pkiv are linear independent points on X. Consider the tangent spaces Ti,...,Tfc^ to 
CX at these points, and let Ui = pT{Ti), i = 1, kw be the corresponding linear spaces in 
the tangent cone to X* at [H], where Pt is the map induced by p on tangent spaces. Then 
the intersection C/i fl • • • fl Uky^r contains the tangent space to W at [H]. But Pi e U:^, so the 
orthogonal complement of the intersection satisfies 

(C/i n • • • n C/fe^)^ = span(C/i^ U ... U C/fe^^) D span(pi,...,pfe^). 
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We conclude that the plane UiH - ■ - H Uk^ has codimension at least kw at [H], and therefore 
the variety W has codimension at least kw in (CP")^. Since kw > k the lemma follows. □ 

At this point, we are pretty much done. To recap, recall that we have shown Y* C X^''\ 
dim(F*) = n — k and dim(X['^]) < n — k. If Xl'^l is irreducible, then we have Y* = Xt'^l and 
Y = (X^^l)* follows. Otherwise, if Xt'^l has several components, then its dual (X''^!)* is the 
union of the dual varieties of each component. One of these components is Y, and hence Y* 
is a component of (X''^!)*. Therefore, the formula fll.ip in Theorem 11.11 is indeed true. □ 

Theorem 11.11 extends in a straightforward manner to reduced and reducible compact real 
algebraic sets with isolated singularities. A colorful picture of a trigonometric space curve 
X with a singularity on the boundary of P = conv(X) is shown in [151 Figure 6]. Also, 
in Subsection 4.1 below we shall examine a reducible space curve with isolated singularities 
with the property that some (finitely many) hyperplanes that are tangent at infinitely many 
points. 

Let X C be a finite union of compact varieties, and assume that X has only isolated 
singularities. As before, we write X be its Zariski closure in CP". For any positive integer 
k we now take X''^! to be the Zariski closure in (CP")^ of the set of all hyperplanes that are 
tangent to X a.t k — s regular points and pass through s singularities on X, for some s, such 
that the (k — s) + s = k points span a (A;— l)-plane. Thus X^^^ contains the dual variety, but, 
in addition, it also contains a hyperplane for each isolated singularity of X. We consider, as 
above, the nested chain of projective varieties 

c ■ ■ ■ c x'2] c xW C (CP")^. 

The algebraic boundary of P = conv(X) is dual to the various X^^l appearing in this chain: 

Theorem 3.2. Let X be a finite union of compact real algebraic varieties that affinely spans 
M", and assume that X has only isolated singularities and that only finitely many hyperplanes 
in CP" are tangent to X at infinitely many points. The algebraic boundary of its convex hull, 
P = conv(X), is computed by biduality using the same formula U.l\) as in Theorem \l.l\ In 
particular, every irreducible component of daP is a component of [X^^^]* for some k. 

Proof. Following the argument of the proof of Theorem II. 1^ we first note that a hyperplane 
H that supports a (A; — l)-dimensional face of P must intersect X in points that span 
the face. Furthermore, H must be tangent to X at the smooth intersection points. Let Y 
be an irreducible component having full-dimensional intersection with the boundary dP of 
conv(X), and let g be a general smooth point on dP H Y . In the notation of the above 
proof, a local parameterization of F at g will involve singular points pi, ...,ps and smooth 
points Ps+i, ■■■,Pk- The /c-tuples (pi, . . . ,pk) of points whose linear span is contained in Y 
form a subvariety Z in the Cartesian product X'^. Since the singular points are isolated, we 
may assume that the restriction of Z to the first s factors is a point. The hypersurface Y is 
therefore a cone with vertex containing the s singular points. The tangent hyperplane to Y 
at q contains the vertex and the tangent spaces at the k — s smooth points, so it coincides 
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with the supporting hyperplane H. The latter part of the proof of Theorem 11.11 applies 
directly to arrive at the same conclusion. □ 

At present, we do not know how to extend our formula (11.11) for the algebraic boundary 
to the convex hull of a compact real variety X whose real singular locus has dimension > 1. 
Also, we do not yet know how to remove the hypothesis that only finitely many hyperplanes 
are tangent to X at infinitely many points. This issue is related to the study of degeneracies 
in [1] and we hope that the techniques introduced in that paper will help for our problem. 

4. More Examples and Computational Thoughts 

We further illustrate our formula for the algebraic boundary of the convex hull of a real 
variety with three concrete examples, starting with a curve that is reducible and singular. 

4.1. Circles and spheres in 3-Space. Let = 3 and suppose that X = Ci U C2 U ■ ■ ■ U C,. 
is the reducible (and possibly singular) curve obtained by taking the union of a collection of 
r > 3 sufficiently general circles Cj that lie in various planes in M^. We have 

(4.1) daP c (X^^l)* u (x[=^])*. 

The surface (Xl^l)* is the union of planes that are tangent to three of the circles and planes 
spanned by the circles. The edge surface (X'^^)* decomposes into quadratic surfaces, namely, 
its components are cylinders formed by stationary bisecant lines defined by pairs of circles. 

For a concrete configuration, consider the convex hull of r = 4 pairwise touching circles 
in 3-space. The surface (X^^l)* is a union of six cylinders, each wrapped around two of the 
circles, while (X^^l)* is the union of planes tangent to three of the circles (four of which contain 
the fourth circle). A picture of this 3-dimensional convex body P is shown in Figure [H Its 
boundary consists of 6 + (4 + 4) = 14 distinct surface patches, corresponding to the pieces in 
(14. ip . which holds with equality. There are six cylinders, four planes touching exactly three 
of the circles, and four planes containing one of the circles and touching the three others. 

A nice symmetric representation of the curve X = Ci U C2 U C3 U C4 is given by the ideal 

{acgt , + + g"^ + f - 2ac-2ag -2at- 2cg - 2ct - 2gt), 

where the variety of that ideal is to be taken inside the probability simplex 

A3 = {(a,c,^,t) eMto : a + c + ^ + t = l}. 

The convex body P looks combinatorially like a 3-polytope with 18 vertices, 36 edges and 20 
cells. Eight of the 20 cells are fiat facets. First, there are the planes of the circles themselves. 
For instance, the facet in the plane t = is the disk {a^ + + (7^ < 2ac + 2ag + 2cg} in the 
triangle {a-\- c-\- g = 1}. Second, there are four triangle facets, formed by the unique planes 
that are tangent to exactly three of the circles. The equations of these facet planes are 

= -a + 2c + 2g + 2t, P^ = 2a-c + 2g + 2t, 
Pg = 2a + 2c-g + 2t, Pt = 2a + 2c + 2g- t. 
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Figure 1. Schlegel diagram of the convex hull of four pairwise touching circles 

The remaining 12 cells in dP are quadratic surface patches that arise from the pairwise 
convex hull of any two of the four circles. This results in 6 quadratic surfaces each of which 
contributes two triangular cells to the boundary. The equations of these six surfaces are 



Qac 


a" 


+ 0^ 


+ 9' 


+ e 


+ 2(ac — ag — eg — at — ct — gt) , 


Qag 


a2 


+ c2 


+ 9' 


+ e 


— 2(ac — ag + eg + at + ct + gt) , 


Qag 


a2 




+ 9' 


+ e 


— 2(ac + ag — eg + at + ct + gt) , 


Qcg 


a2 


+ c2 


+ 9' 


+ e 


— 2(ac + ag + eg — at + ct + gt) , 


Qct = 


a2 


+ c2 


+ 9' 


+ e 


— 2(ac + ag + eg + at — ct + gt) , 


Qgt = 


a" 


+ c2 


+ 9' 


+ e 


— 2(ac + ag + eg + at + ct — gt) . 



Each circle is subdivided into six arcs of equal length. Three of the nodes arise from inter- 
sections with other circles, and the others are the intersections with the planes Pa-, Pc, Pg,Pt- 
This accounts for all 18 vertices and 24 "edges" that are arcs. The other 12 edges of dP 
are true edges: they arise from the four triangles. These are shown in green in the Schlegel 
diagram of Figured] The 12 cells corresponding to the six quadratic surfaces are the 12 ruled 
cells in the diagram, and they come in pairs according to the six different colors. The six 
intersection points among the 4 circles are indicated by black dots, whereas the remaining 
twelve vertices correspond to the green dots which are vertices of our four green triangles. 
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Figure 2. The Zitms surface + 2;^ + (y^ - 1)^ = 

4.2. Zitrus. We have seen that the convex hull of algebraic surfaces in can have surfaces 
of very high degree in its boundary. For instance, if X is a general smooth surface of degree 
d = Q then the bitangent surface (Xt^')* has degree 3168. On the other hand, that number 
can be expected to drop substantially for most singular surfaces. Let us consider the sextic 

f{x,y,z) = x' + z^ + {y'-l)^ 

The surface X = V{f) in is taken from Herwig Hauser's beautiful Gallery of Algebraic 
Surfaces. The name given to that surface is Zitrus. It appears on page 42-43 of the catalog [9J 
of the exhibition Imaginary. For an electronic version see www . f reigeist . cc/gallery . html. 

We choose affine coordinates (a, b, c) on the space of planes ax + by + cz + 1 = inM.^. In 
these coordinates, the variety X^^l is the union of two quadratic curves given by the ideal 

{b + 1, 27a^ + 27c^ - 16) n (6-1, 27a^ + 27c^ - 16). 

These curves parametrize the tangent planes that pass through one of the two singular points 
of the Zitrus. Each curve dualizes to a singular quadratic surface, and (X'^^)* is given by 

(16x^ - 27y^ + 16^2 + 54?/ - 27) U {I6x^ - 27y^ + 16z'^ - 5Ay - 27). 

The Zitrus X has no tritangent planes, so daP = X U (X^^l)*, and we conclude that the 
algebraic boundary of the convexified Zitrus P = conv(X) has degree 10 = 6 + 2 + 2. 
We now perturb the polynomial / and consider the smooth surface X = V{f) defined by 

f{x,y,z) = x^ + z^ + {y^-lf-l. 
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The curve of bitangent planes, X^'^\ has again two components. It is defined by the ideal 
( 6 , + c2 - 1 ) n ( 900^62 _ gefo^ + 90^2^2 _ i29a2 + 12862 - 129c2 + 48 , 

135a^ - 1446^ + 2700^02 + 135c^ - 6a^ + 272b^ - Gc^ - 48 ). 

The first curve dualizes to the cylinder {x"^ + 2"^ = !}• The other component of the boundary 
surface (X'21)* has degree 16. Its defining polynomial has 165 terms which start as follows: 

16777216x^^-169869312x^^?/2 + 1601372160x^V-7081205760a;^V + 26435102976xV 

4.3. Grassmannian. We consider the oriented Grassmannian X = Gr(2, 5) of oriented two- 
dimensional linear subspaces of MP . This is the 6-dimensional subvariety of defined by 

( P12 + Pl3 + Pl4 + Pl5 + P23 + pIa + P25 + pIa + P35 + - 1' P12P34 - P13P24 + P14P23, 
P12P35-P13P25+P15P23, P12PA5-PiaP25+P15P2A, Pl3PA5-PlAP35+Pl5P3A, P23PA5-P2AP35+P25P3a) ■ 

Its convex hull P = conv(X) is a Grassmann orbitope, a class of convex bodies that are 
of interest to differential geometers. We refer to [12], [T71 §7], and the references given 
therein. The determinant of the Hermitian matrix in the spectrahedral representation of 
P in [T71 Theorem 7.3] has degree 8 and it factors into two quartic factors. Only one of 
these two factors is relevant for us, and we display it below. Namely, the algebraic boundary 
daP = (Xt^l)* is the irreducible hypersurface of degree 4 represented by the polynomial 

Pt2 + Pt3 + P\a + P\^ + P\3 + P\a + P\^ + pIa + Pk + P45 
+2P?2P?3 + 2p?2P?4 + 2p?3P?4 + 2p?2P?5 + '^pWi^ + '^pIaPI^ + M2PI3 + 2p?3P23 " 2PuP23 
-2P?5P23 + 2P?2P24 - 2P?3P24 + '^pIaPIa ' M^pIa + 2^23^24 + 2p?2P25 ' '^P\3PI^ ' MaPI^ 

+'^p1^pI^ + 2pLpi5 + 2P24P25 - 2p?2P34 + 2p?3Pi4 + '^pIaPIa - '^pIA + 2p^3Pi4 + '^pIaPIa 

-2^25^14 - 2p?2Pi5 + 2p?3P35 " '^pIaPI^ + '^pI^pI^ + 2p^3P35 " 2^14^35 + 2^15^35 + 2^14^15 

-2p?2Pl5 - 2p?3Pl5 + MaPI^ + 2p?5Pl5 - 2P23PI5 + '^pIaPI^ + '^pWa^ + 2pLp15 + 2P35PI5 

+8pi3Pl4P23P24 + 8P13P15P23P25 + 8pi4Pl5P24P25 " 8^12^14^23^34 + 8^12^13^24^34 
-8P12P15P23P35 + 8P12P13P25P35 + 8pi4Pl5P34P35 + 8P2AP25P3AP35 " 8P12P15P2APA5 
+8P12P14P25P45 - 8pi3Pl5P34P45 - 8P23P25P3APA5 + 8^13^14^35^45 + 8^23^24^35^45 

-2^2^ - 2^23 - 2pl - 2pl, - 2pl, - 2pl, - 2pl, - 2pl, - 2pl, - 2pl, + 1. 

This hypersurface represents a 6-dimensional family of 3-dimensional facets of P. Each facet 
of P is a 3-dimensional ball. It meets the variety X in its boundary, which is a 2-sphere. 

4.4. Computing. This paper raises the following algorithmic problem: given a projective 
variety X, either by its ideal or by a parametrization, how to compute the equations defining 
[X^^^y in practise? The passage from X to X^^^ can be phrased as an elimination problem in a 
fairly straightforward manner. In principle, we can use any Grobner-based computer algebra 
system to perform that elimination task. However, in our experience, this approach only 
succeeds for tiny low-degree instances. Examples such as the Grassmannian in Subsection 
14.31 appear to be out of reach for a general purpose implementations of our formula (11. ip . 

Even the first instance = 1, which is the passage from a variety X to its dual variety 
X*, poses a considerable challenge for current computational algebraic geometry software. 
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The case of plane curve is still relatively easy, and it has been addressed in the literature [5] . 
However, what we need here is the case when X is not a hypersurface but X* is. The first 
interesting situation is that of a space curve X C CP^. Our computations for space curves, 
both here and in [13], were performed in Macaulay2 [8], but, even with ad hoc tricks, they 
turned out to be more difficult than we had expected when we first embarked on our project. 
Here is an illustration of the issue. Let X be the smooth sextic curve in defined by 

/ 2 I 2 I 2 I 2 3 \ 

\x + y +z +w,xyz — w ). 

The following lines of Macaulay2 code find the surface X* in (CP^)^ that is dual to X: 
S = QQ[x,y,z,w,X,Y,Z,W] ; 

d = 4; pairing = first sum(d, i->(gens S)_i*(gens S)_-[i+d}-); 
makedual = I -> (e = codim I ; J = 

saturated + minors (e+l,submatrix(jacobian(I+ideal (pairing) ) ,-[0. .d-l},)), 
minors (e , submatrix( jacobian(I) ,-[0 . . d-l} , ) ) ) ; eliminate ( (gens S)_-[0 . .d-l} , J) ) 
makedual ideal ( x~2+y~2+z"2+w"2 , x*y*z-w"3 ); 

This program runs for a few minutes and outputs a polynomial of degree 18 with 318 terms: 

729a;^^?/^ + 3861x^'^y^ + 7954a;^°?/® + 7954a;^?/i° + 3861x%^^ + 729xY'^ + UbSx^^z"^ + ■■■ 

Projective duality tends to produce large equations, even on modestly sized input, and 
symbolic programs, like our little Macaulay2 fragment above, will often fail to terminate. 

One promising alternative line of attack is offered by numerical algebraic geometry [2]. 
Preliminary experiments by Jonathan Hauenstein demonstrate that the software Bertini 
can perform the transformations X X* and X H- (X''^')* in a purely numerical manner. 

Convex algebraic geometry requires the development of new specialized software tools, 
both symbolic and numeric, and integrated with optimization method. The advent of such 
new tools will make our formula (11. ip more practical for non-linear convex hull computations. 
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